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The late infall of cold dark matter onto an isolated galaxy such as our own produces 
flows with definite velocity vectors at any physical point in the galactic halo. It also 
produces caustics which are places where the dark matter density is very large. The 
outer caustics are topological spheres whereas the inner caustics are rings. The 
self-similar model of galactic halo formation predicts that the caustic ring radii 
a n follow the approximate law a n ~ 1/ra. In a recent study of 32 extended and 
well-measured galactic rotation curves, we found evidence for this law. 



1 Introduction 

Before the onset of galaxy formation but after the time t eq of equality between 
matter and radiation, the velocity dispersion of the cold dark matter candidates 

is very small, of order 5v a (t) ~ 3 • 10~ 17 ( 10 m eV i}f^} 2 ^ f° r axions and 

1/2 

Sv w (t) ~ 1CT 11 (¥) 2/3 for WIMPs, where t is the present age of 

the universe and m a and my/ are respectively the masses of the axion and 
the WIMP. In the context of galaxy formation, such small velocity dispersions 
are entirely negligible. Massive neutrinos, on the other hand, have primordial 
velocity dispersion dv u (t) ~ 5.3 10~ 4 (fpf^) (y) 2 ^ 3 which is comparable to 
the virial velocity- in galaxies and therefore non-negligible in the context of 
galaxy formation El. This is the reason why massive neutrinos are called 'hot 
dark matter'. 

Collisionless dark matter particles lie on a thin 3-dimensional (3D) sheet 
in 6D phase-space. The thickness of this sheet is the primordial velocity dis- 
persion Sv. If each of the aforementioned species of collisionless particles is 
present, the phase-space sheet has three layers, a very thin layer of axions, a 
medium layer of WIMPs and a thick layer of neutrinos. The phase-space sheet 
is located on the 3D hypersurface of points (r, v) : v = H(t)r + Av(r, t) where 
H{t) = is the Hubble expansion rate and Av(r, t) is the peculiar velocity 
field. Fig. 1 shows a 2D section of 6D phase-space along the (z, z) plane. 
The wiggly line is the intersection of the 3D sheet on which the particles lie 
in phase-space with the plane of the figure. The thickness of the line is the 
velocity dispersion Sv, whereas the amplitude of the wiggles in the line is the 
peculiar velocity Av. If there were no peculiar velocities, the line would be 
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Figure 1: The wiggly line represents the intersection of the (z, z) plane with the 3D sheet 
on which the collisionless dark matter particles lie in phase-space. The thickness of the line 
is the primordial velocity dispersion. The amplitude of the wiggles in the i direction is the 
velocity dispersion associated with density perturbations. Where an overdensity grows in 
the non-linear regime, the line winds up in clockwise fashion. One such overdensity is shown. 



straight since z — H(t)z in that case. 

The peculiar velocities are associated with density perturbations and grow 
by gravitational instability as Av ~ i 2 / 3 . On the other hand the primor- 
dial velocity dispersion decreases on average as Sv ~ i~ 2 / 3 , consistently with 
Liouville's theorem. When an overdensity enters the non-linear regime, the 
particles in its vicinity fall back onto it. This implies that the phase-space 
sheet 'winds up' there in clockwise fashion. One such overdensity is shown in 
Fig. 1. In the linear regime, there is only one value of velocity, i.e. one single 
flow, at a typical location in physical space, because the phase-space sheet cov- 
ers physical space only once. On the other hand, inside an overdensity in the 
non-linear regime, the phase-space sheet covers physical space multiple times 
implying that there are several (but always an odd number of) flows at such 
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locations. 

At the boundary surface between two regions one of which has n flows 
and the other n + 2 flows, the physical space density is very large because 
the phase-space sheet has a fold there. At the fold, the phase-space sheet is 
tangent to velocity space and hence, in the limit of zero velocity dispersion 
(Sv = 0), the density diverges since it is the integral of the phase-space density 
over velocity space. The structure associated with such a phase-space fold is 
called a 'caustic'. It is a surface in physical space. It is easy to show that, in 
the limit of zero velocity dispersion, the density diverges as d ~ when the 
caustic is approached from the side with n + 2 flows, where a is the distance 
to the caustic. Velocity dispersion cuts off the divergence. 

As mentioned above, the process of galactic halo formation involves the 
local winding up of the phase-space sheet of collisionless dark matter particles. 
If the galactic center is approached from an arbitrary direction at a given time, 
the local number of flows increases. First, there is one flow, then three flows, 
then five, seven... The number jpf flows at our location in the Milky Way 
galaxy today has been estimated □ to be of order 100. The boundary between 
the region with one (three, five, ...) and the region with three (five, seven, ...) 
flows is the location of a caustic which is topologically a sphere surrounding the 
galaxy. When these caustic spheres are approached from the inside the density 
diverges as d ~ in the zero velocity dispersion limit. These spheres are the 
outer caustics in the phase-space structure of galactic halos. In addition there 
are inner caustics. 

It is a little more difficult to see why there must be inner caustics, and to 
derive their structure. See refs. [3,4] for details. The inner caustics are rings. 
They are located near where the particles with the most angular momentum in 
a given in and out flow are at their distance of closest approach to the galactic 
center. A ring is a closed tube whose cross-section is a £>_4 catastopheQ. The 
cross-section is shown in Fig. 2 in the limit of axial and reflection symmetry, 
and where the transverse dimensions, p and q, are much smaller than the ring 
radius a. In the absence of any symmetry, the cross-section of the tube does not 
have the exact shape shown in Fig. 2 but it still has that shape qualitatively, 
i.e. it is still a closed line with three cusps one of which points away from the 
galactic center. 

The existence of caustic rings of dark matter follows from only two as- 
sumptions: 

1. the existence of collisionless dark matter 

2. that the velocity dispersion of the infalling dark matter is much less, by 
a factor ten say, than the rotation velocity of the galaxy. 



3 



p 



Figure 2: Cross-section of a caustic ring in the case of axial and reflection symmetry. The 
galactic center is to the left of the figure. The z-dircction is orthogonal to the galactic plane. 
The p-direction is radial, a is the caustic ring radius. The closed line with three cusps is the 
location of a caustic surface. The density diverges when the surface is approached from the 
inside as cr -1 / 2 where a is the distance to the surface. 



Only the second assumption requires elaboration. Velocity dispersion has the 
effect of smoothing out caustics. The question is when is the velocity disper- 
sion so large as to smooth caustic rings over distance scales of order the ring 
radius a, thus making the notion of caustic ring meaningless. In ref. [3] this 
critical velocity dispersion was estimated to be 30 km/s = 10~ 4 for the caustic 
rings in our own galaxy, whose rotation velocity is 220 km/s. 10~ 4 is much less 
than the primordial velocity dispersion Sv of the cold dark matter candidates. 
However the velocity dispersion Aw associated with density perturbations also 
smoothes caustics in coarse grained observations. So the question is whether 
the velocity dispersion Av of cold dark matter particles associated with den- 
sity perturbations falling onto our galaxy is less than 30 km/s. The answer 
is yes with high probability since the infalling dark matter particles are not 
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associated with any observed inhomogeneities. 30 km/s is of order the velocity 
dispersion of the Magellanic Clouds. For the velocity dispersion of the dark 
matter particles presently falling onto our galaxy to be as large 30 km/s, these 
particles would have to be part of clumps whose mass/size ratio is of order that 
of the Magellanic Clouds. But if that were the case, why did these clumps fail 
to become luminous? 

One might ask whether caustic rings can be seen in N-body simulations of 
galaxy formation. The generic surface-caustics associated with simple folds of 
the phase-space sheet have been seenB. However, caustic rings would require 
far greater resolution than presently available, at least in a 3D simulatioii of 
our own halo. Indeed, the largest ring in our galaxy has been estimated^ to 
have radius of order 40 kpc. It is part of an in and out flow that extends to 
the Galaxy's current turnaround radius, of order 2 Mpc. To resolve this first 
ring, the spatial resolution would have to be considerably smaller than 10 kpc. 
Hence a minimum of 2 • n^^ys -j- (2Mpc) 3 ~ 7 • 10 7 particles would be required 
to see the caustic ring in a simulation of this one flow. However, the number 
of flows at 40 kpc in our halotJ is of order 10. So it appears that 10 9 particles 
is a strict minimum in a 3D simulation of our halo. Even so, this addresses 
only the kinematic requirement of resolving the halo in phase-space, assuming 
moreover that the particles are approximately uniformly distributed on the 
phase-space sheets. There is a further dynamical requirement that 2-body 
collisions do not artificially 'fuzz up' the phase-space sheets. Indeed 2-body 
collisions are entirely negligible in the flow of cold dark matter particles such 
as axions or WIMPs. On the other hand, 2-body collisions are present, and 
hence the velocity dispersion is artificially increased, in the simulations. This 
may occur to such an extent that the caustics are washed away even if 10 9 
particles are used. 

In the self-similar infall model Qa of galactic halo formation the caustic 
ring radii a n are predicted 0: 

W : „ = 1,2,..) a (39, 19.5, 13. 10. 8, ..Jkp. (gg) (H) 

(1) 

where h is the present Hubble rate in units of 100 km/s. Mpc, v rot is the rotation 
velocity of the galaxy and j ma x is the maximum of its dimensionless angular 
momenturcLjdistribution as defined in ref. [8]. In Eq. (1) we assume that the 
parameterQa e = 0.3. 

Eq. (Q) predicts the caustic ring radii of a galaxy in terms of its first ring 
radius a\. If the caustic rings lie close to the galactic plane they cause bumps 
in the rotation curve, at the caustic ring radii. As a possible example of this 



5 



effect, consider^ the rotation curve of NGC3198, one of the best measured. It 
has three faint bumps at radii: 28, 13.5 and 9 kpc, assuming h — 0.75. The 
ratios happen to be consistent with Eq. ([!]) assuming the bumps are caused 
by the first three (n — 1,2,3) ring caustics of NGC3198. Moreover, since 
frot = 150 km/s, jmax is determined to equal 0.28. The uncertainty in h is 
a systematic effect that can be corrected for when determining j max because 
the bump radii scale like l/h' where h! is the Hubble rate assumed by the 
observer in constructing the rotation curve, and the caustic ring radii scale as 
/h. Rises in the inner rotation curve of the Milky Way were also interpreted 
as due to caustics n = 6, 7, 8, 9, 10, 11, 12 and 13. This determined the value 
of jmax of our own galaxy to be 0.263. The first five caustic ring radii in our 
galaxy are then predicted to be: 41, 20, 13.3, 10, 8 kpc. 



2 Evidence for universal structure in galactic halos 

Motivated by these findings, we analyzed! a set of 32 extended well-measured 
galactic rotation curves which had been previously selectedE3 under the criteria 
that each is an accurate tracer of the galactic radial force law, and that it 
extends far beyond the edge of the luminous disk. 

According to the self-similar caustic ring model, each galaxy has its own 
value of jmax- Over the set of 32 galaxies selected in ref. [10], j m ax has some 
unknown distribution. However, the fact that the values of j ma x of NGC3198 
and of the Milky Way happen to be close to one another, within 7%, suggests 
that the jmax distribution may be peaked near a value of 0.27 . Our strategy 
is to rescale each rotation curve according to 

. 220 km/s, 

r^f = r{ '—) (2) 

Vrot 

and to add them in some way. Since Eq. (Q) predicts the nth caustic radius 
a n to be distributed like j ma x for all n, and it fixes the ratios a n /a± ~ the 
sum of rotation curves should show the j m ax distribution, once for n = 1, then 
at about half the n = 1 radii for n — 2, then at about 1/3 the n = 1 radii for 
n = 3, and so on. If the jmax distribution is broad, the sum of rotation curves 
is unlikely to show any feature. However, if it is peaked, then the sum should 
show a peak for n — 1 at some radius, then again at 1/2 that radius for n = 2, 
at 1/3 the radius for n — 3, and so on. If the j ma x distribution is peaked at 
0.263 (the value for the Milky Way) the peaks in the sum of rotation curves 
should appear at 41 kpc, 20 kpc, 13.3 kpc ... 

The procedure followed to add the 32 rotation curves is described in detail 
in ref. [9]. Briefly, we proceeded as follows. For each rotation curve, all data 
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Figure 3: Binned data for 32 galaxy sample, with peaks fit to Gaussians. 

points with rescaled radii f < 10 kpc were deleted to remove the effect of the 
luminous disk. The remaining points were then fitted to a line. The rotation 
velocity v ro t used to rescale the radii in Eq. (0) is the average of that line. The 
rms deviation y/ (Sv 2 ) from the linear fit was determined for each galaxy. This 
was taken to be the error on the residuals Svi, i.e. the differences between the 
measured velocities in a rotation curve and the linear fit. Finally the sample 
of 32 galaxies was averaged in 2 kpc radial bins: 

1 Nl 

1 J=l 

where N% is the number of data points in the bin. The assigned error on each 
hi is then simply \]\fN\. Fig. 3 shows the result. 

There are two features evident at roughly 20 and 40 kpc. A fit to two 
Gaussians plus a constant indicates features at 19.4±0.7kpc and 41.3±0.8kpc, 
with overall significance of 2.4ct and 2.6er, respectively. Fig. 3 shows the fitted 
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curve. When the same fit is applied to the same data in 1 kpc bins, the 
significance of the two peaks is 2.6 and 3.0 a respectively. The locations of 
the features agrees with the predictions of the self-similar caustic ring model 
with the jmax distribution peaked at 0.27. The use of Gaussians to fit the 
peaks in the combined rotation curve was an arbitrary choice in the absence 
of information on the j max distribution. 

The existence of velocity peaks and caustic rings in_the cold dark matter 
distribution is relevant to axiontil and WIMP searchpdlj. Caustics may also 
be investigated using gravitational lensing techniques c3. 
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